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The muon capture reactions 2 H(/i~, v^)nn and 3 He(/i~ , z/ M ) 3 H are studied with conventional or 
chiral realistic potentials and consistent weak currents. The initial and final A = 2 and 3 nuclear 
wave functions are obtained from the Argonne Uig or chiral N3LO two-nucleon potential, in com- 
bination with, respectively, the Urbana IX or chiral N2LO three-nucleon potential in the case of 
A = 3. The weak current consists of polar- and axial-vector components. The former are related 
to the isovector piece of the electromagnetic current via the conserved-vector-current hypothesis. 
These and the axial currents are derived either in a meson-exchange or in a chiral effective field 
theory (^EFT) framework. There is one parameter (either the iV-to-A axial coupling constant in 
the meson-exchange model, or the strength of a contact term in the ^EFT model) which is fixed by 
reproducing the Gamow- Teller matrix element in tritium /3-decay. The model dependence relative 
to the adopted interactions and currents (and cutoff sensitivity in the xEFT currents) is weak, 
resulting in total rates of 392.0 ± 2.3 s _1 for A = 2, and 1484 ± 13 s _1 for A = 3, where the spread 
accounts for this model dependence. 

PACS numbers: 23.40.-s,21.45.-v,27.10.+h 



I. INTRODUCTION 

There is a significant body of experimental and theo- 
retical work on muon capture in nuclei (see Refs. [l|, 0] 
for a review). These processes provide a testing ground 
for wave functions and, indirectly, the interactions from 
which these are obtained, and for models of the nuclear 
weak current. This is particularly important for neu- 
trino reactions in light nuclei [3] and for processes, such 
as the astrophysically relevant weak captures on proton 
and 3 He, whose rates cannot be measured experimen- 
tally, and for which one has to rely exclusively on theory. 
Thus, it becomes crucial to study within the same theo- 
retical framework related electroweak transitions, whose 
rates are known experimentally Muon captures are 
among such reactions. 

In the present work, we focus our attention on muon 
capture on deuteron and 3 He, i.e., on the reactions 



3 He ->■ 3 H 



(1.1) 
(1.2) 



Muon capture on 3 He can also occur through the two- 
(nd) and three-body (nnp) breakup channels of 3 H. How- 
ever, the branching ratios of these two processes are 20% 
and 10%, respectively, and will not be considered in the 
present work. 

These reactions have been studied extensively through 
the years, experimentally and theoretically. In reac- 
tion p.l[) . the stopped muons are captured from two hy- 
pcrfine states, / = 1/2 or 3/2. The doublet capture rate 
r D has been calculated by severalgroups to be about 40 
times larger than the quadruplet [lj [|J . We will therefore 
consider only T D . The first attempt to measure T D was 



carried out over forty years ago by Wang et al.: they ob- 
tained T D = 365(96) s" 1 [|. A few years later, Bertin 
et al. measured T D — 445(60) s" 1 |6j. Measurements 
performed in the eighties gave T D = 470(29) s" 1 0] and 
T D = 409(40) s _1 8]. These measurements, while con- 
sistent with each other, are not very precise - errors are 
in the 6-10 % range. However, there is hope to have 
this situation clarified by the MuSun Collaboration Q, 
which is expected to perform an experiment at the Paul 
Scherrer Institut, with the goal of measuring T D with a 
precision of 1 %. 

The experimental situation for reaction (|1.2|) is much 
clearer: after a first set of measurements in the early 
sixties a very precise determination in the late 

nineties yielded a total capture rate r =1496(4) s _1 
a value consistent with those of the earlier measurements, 
although these were affected by considerably larger un- 
certainties. 

Theoretical work on reactions (jl.l[) and (|1 .2[) is just as 
extensive, and a list of publications, updated to the late 
nineties, is given in Table 4.1 of Ref. [lj], and in Ref. Q. 
Here, we limit our considerations to the calculations of 
Refs. (l5T - fl7l|. W e also comment on the recent studies of 
Ando et al. 18| and Ricci et al. [h| . 

The calculations of Refs. [l5T - [l~7T | were performed 
within the "Standard Nuclear Physics Approach" 
(SNPA): their authors used the realistic potential mod- 
els available at the time to obtain the nuclear wave func- 
tions, and included in the nuclear weak current operator 
one-body (impulse approximation) and two-body opera- 
tors. In Ref. [13, T D was calculated to be 416(7) s" 1 , the 
uncertainty coming from imprecise knowledge of the cou- 
pling constants, and in Refs. [l6[ and fl7l 3 99 s _1 and 402 
s _1 , respectively. The results of Refs. |l6| and [l7j are in 
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good agreement with each other, while that of Ref. [l5[ 
differs by ~ 4 %. It is important to note, however, that 
the meson-exchange currents (MEC) contributions were 
not constrained to reproduce any experimental observ- 
able, such as the triton half-life, as is now common prac- 
tice din. 

Reference [l8[ overcame some of the limitations in- 
herent to the earlier studies. Much along the lines of 
Ref. [22j , reaction (ll.lj) was studied within a hybrid chiral 
effective field theory (^EFT) approach, in which matrix 
elements of weak operators derived in %EFT were eval- 
uated between wave functions obtained from a realistic 
potential, specifically the Argonne Vi$ (AV18) [23|. The 
^EFT axial current contains a low-energy constant which 
was fixed by reproducing the experimental Gamow- Teller 
matrix element (GT EXP ) in tritium /3-decay. The cal- 
culation, however, retained only the S*-wave contribu- 
tion in the nn final scattering state (the 1 Sq state), 
and higher partial-wave contributions were taken from 
Ref. [l6|. This approach yielded a value for T D of 386 
s" 1 , with r D ( 1 S u )=245(l) s^ 1 , the theoretical error be- 
ing related to the experimental uncertainty in GT EXP . 

The latest SNPA calculation of muon capture on 
deuteron has been carried out in Ref. and has led 
to values in the range of 416-430 s _1 (see Table 1 of 
Ref. depending on the potential used, the Nijmegen 

I or Nijmegen 93 [24j. However, the model for the ax- 
ial current is not constrained by data, resulting in the 
relatively large spread in T D values. 

Finally, there is a calculation based on pionless 
EFT [25] with the objective of constraining the two- 
nucleon axial current matrix element by reproducing the 
muon capture rate on deuterons. This same matrix ele- 
ment enters the pp weak capture. 

Theoretical studies for reaction (ll.2p within the SNPA 
have been performed in the early nineties by Congleton 
and Fearing [26| and Congleton and Truhlik (2?J ■ In this 
later work, the nuclear wave functions were obtained from 
a realistic Hamiltonian based on the Argonne U14 (AV14) 
two-nucleon [28| and the Tucson-Melbourne (TM) three- 
nucleon [29j interactions. The nuclear weak current re- 
tained contributions similar to those of Ref. [l9j . The 
value obtained for the total capture rate To was 1502(32) 
s _1 , the uncertainty due to poor knowledge of some of 
the coupling constants and cutoff parameters entering the 
axial current. 

A first study of the model-dependence of predictions 
for the total rate of muon capture on 3 He was carried out 
in Ref. 30] , within the SNPA, but without the inclusion 
of MEC contributions. Values for To were found to vary 
by ~ 100 s _1 , depending on the potential model consid- 
ered. Both old- (Paris I31J, Bonn A and B [32|, AV14) 
and new-generation (Nijmegen I [24j and CD-Bonn [331 ]) 
potentials were used. However, no three-nucleon forces 
were included. This second aspect of the calculation, 
along with the absence of MEC contributions, might be 
the origin of the large model dependence observed in the 
results. 



A first attempt to study muon capture on 3 He in 
a way that was consistent with the approach adopted 
for the weak proton capture reactions, involved some of 
the present authors [34j . The nuclear wave functions 
were obtained with the hyperspherical-harmonics (HH) 
method (see Ref. [35| for a recent review) , from a realistic 
Hamiltonian based on the AV18 two-nucleon and Urbana 
IX [36| (UIX) three-nucleon interactions. The model for 
the nuclear weak current was taken from Refs. [2(| [2l| . 
However, two additional contributions were included: 
the single-nucleon pseudoscalar charge operator and the 
pseudoscalar two-body term in the iV-to-A transition ax- 
ial current. Both contributions are of order 0(q 2 /m 2 ), 
where q is the momentum transfer in the process and m 
is the nucleon mass, and were obviously neglected in the 
pp and hep captures of Refs. [20l l2lj]. for which q << m. 
The axial coupling constant for the iV-to-A transition 
was constrained to reproduce GT EXP . The total cap- 
ture rate To was found to be 1484(8) s^ 1 , where the un- 
certainty results from the adopted fitting procedure and 
experimental error on GT EXP . A calculation based on 
the older AV14/TM Hamiltonian model yielded a To of 
1486(8) s _1 , suggesting a weak model-dependence. 

Recently, a hybrid calculation has appeared [37|, m 
which the nuclear wave functions have been obtained 
with the Effective Interaction HH method [38[, and the 
%EFT weak current is that of Ref. 22]. It has yielded 
a value for Tq of 1499(16) s _1 , where the error has two 
main sources: the experimental uncertainty on the triton 
half-life, and the calculation of radiative corrections. 

In light of the previous considerations, it is clear that a 
calculation is still lacking which: (i) treats reactions (|1.1[) 
and (|1.2[) simultaneously in a consistent framework, ei- 
ther SNPA or %EFT; (ii) is based on up-to-date Hamilto- 
nian models to generate the wave functions; (iii) reduces 
the model dependence of the weak axial current by con- 
straining it to reproduce GT EXP . The goal of the present 
work is to fill this gap. The calculation has been struc- 
tured as follows: the nuclear Hamiltonian models consid- 
ered consist of the AV18 and N3LO [39] two-nucleon in- 
teractions for the A — 2 systems, augmented by the UIX 
and N2LO 40] three-nucleon interactions for the A = 3 
systems. The nuclear weak current is derived from ei- 
ther SNPA or xEFT. In both cases, its axial component 
is calibrated by fitting GT EXP , while its vector part is 
related to the isovector electromagnetic current by the 
conserved-vector-current (CVC) hypothesis. The SNPA 
version of it reproduces well static magnetic properties 
of few-nucleons systems. In the xEFT electromagnetic 
current of Ref. [42| , adopted in the present work, the two 
low-energy constants are fixed by reproducing the trinu- 
cleon magnetic moments. 

In closing, we note that the study of Ref. [34| has 
established that the total rate of reaction (ll.2p scales 
approximately linearly with the trinucleon binding en- 
ergy. Thus a realistic calculation of this rate must, at a 
minimum, include a three-nucleon potential which repro- 
duces well these energies (as is the case here). Finally, 
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we could have adopted the set of two-nucleon chiral po- 
tentials at N3LO derived by the Bonn group [U] for a 
range of cutoff parameters A, in order to explore the 
sensitivity of the results to the short-range behavior of 
the potentials. However, the potentials considered in the 
present study (AV18 and N3LO) have so drastically dif- 
ferent treatments of this short-range behavior that they 
should provide a meaningful measure of the model de- 
pendence originating from it. 

The paper is organized as follows. In Sec. [TT] we list 
the explicit expressions for the observables of interest in 
terms of reduced matrix elements of multipole operators. 
In the case of reaction (|1.2[) the derivation is given in 
Ref. [34j]. In Sec. IIII| we briefly review the method used 
to calculate the nuclear wave functions and summarize 
the main results for the observables of the nuclear sys- 
tems involved in the reactions of interest. In Sec. IIV1 we 
describe the model for the nuclear weak current, both its 
SNPA and %EFT versions. In Sec. [V] we present and 
discuss the results obtained for the total capture rates 
of reactions (jl.lj) and (|1.2I) . and finally in Sec. IVII we 
summarize our conclusions. 



II. OBSERVABLES 

The muon capture on deuteron and 3 He is induced by 
the weak interaction Hamiltonian [43[ 

Hw = ^J dxZ CT (x)r(x) , (2.1) 

where Gv is the Fermi coupling constant, Gy=l. 14939 
xlO -5 GeV -2 as obtained from an analysis of + to + 

I 



for the muon capture on deuteron, where p is the nn 
I 



/^-decays [44j> and l a and j a are the leptonic and hadronic 
current densities, respectively. The former is given by 

Z ff (x) - e- ik - x u(k t/ ,^) 7CT (I~ 75 )^(x, S/i ) , (2.2) 

where Vv( x ' s a») * s ^ ne ground-state wave function of the 
muon in the Coulomb field of the nucleus in the initial 
state, and w(k^, h v ) is the spinor of a muon neutrino with 
momentum k„, energy E v (=k v ), and helicity h v . While 
in principle the relativistic solution of the Dirac equation 
could be used, in practice it suffices to approximate 



K -> , (2.3) 

since the muon velocity u„ ~ Za < 1 (a is the fine- 
structure constant and Z=l or 2 for deuteron or 3 He, 
respectively). Here ^\ s {x) is the Is solution of the 
Schrodinger equation and, since the muon is essentially 
at rest, it is justified to replace the two-component spin 
state x(sp) with the four-component spinor u(k M , s^) in 
the limit k M — > 0. This will allow us to use standard 
techniques to carry out the spin sum over s M at a later 
stage. 

In order to account for the hyperfine structure in the 
initial system, the muon and deuteron or 3 He spins are 
coupled to states with total spin /, equal to 1/2 or 3/2 
in the deuteron case, and to or I in the 3 He case. The 
transition amplitude can then be conveniently written as 



(2.4) 

I 

relative momentum, and [34j 



T w {f,fz;s 1 ,s 2 ,h u ) = (nn,s 1 ,s 2 ;v,hv\Hw\{lJt,d);f,fe) 

- %^lE^ s ^ ls dl// Z )^(^>^)(* PiS1S2 (nn)|r(q)|* d ( Sd )) , 



T w (fJ z ;s' 3 ,h v ) = ( 3 H,a' 3 ;u,h v \H w \( Ji , 3 Be);f,f t ) 

~ ^™^(LL 3 \ff z )l a ( hu , SM )(*3 H (4)ir(q)l*3He( S 3)) , (2.5) 



for muon capture on 3 He. In Eqs. (|2.4p and (|2.5[) we have defined 

la(K, s M ) = u(k v ,h v )j a (1 -j^ufo^s^) , (2.6) 
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and the Fourier transform of the nuclear weak current 
has been introduced as 



j a (q) 



dxe it i x r(x)^(p(q)J(q)) , 



(2.7) 



with the leptonic momentum transfer q defined as q = 
— k„ ~ — k„. The function ipi s (x) has been factored 
out from the matrix element of j CT (q) between the initial 
and final states. For muon capture on deuteron, %l>f s is 
approximated as [43[ 



iv^r = i^(o)i 2 = 



(a Mm) 3 



(2.8) 



where V'is(O) denotes the Bohr wave function for a point 
charge e evaluated at the origin, and [i^d is the reduced 
mass of the (//, d) system. For muon capture on 3 He, t/>*, 
is approximated as 34] 



IV* 



av|2 
Is I 



= n (2a^3 He ) 3 



(2.9) 



where in this case is the reduced mass of the 

(/i, 3 He) system, and the factor 1Z approximately ac- 



counts for the finite extent of the nuclear charge distri- 
bution [43]. This factor is defined as 



n 



av 12 



(2.10) 



with 



J dxe lcl ' x ipi s (x)p(x) 
Jrfxc^pfj;) 



(2.11) 



where p(x) is the 3 He charge density. It has been calcu- 
lated explicitly by using the charge densities correspond- 
ing to the two Hamiltonian models considered in the 
present study (AV18/UIX and N3LO/N2LO), and has 
been found to be, for both models, within a percent of 
0.98, the value commonly adopted in the literature (43| . 

In the case of muon capture on deuteron, the final state 
wave function is expanded in partial waves as 



, (nn) 



E 

LJJ, 



i L F L * L , (p) (SS Z , LL Z \JJ Z ) <f JJ * (p) 



(2.12) 



r 



where ^ nn " (j>) is the nn wave function - it will be dis- 
cussed in Sec. IIIII In the present work, we restrict our 
calculation to J < 2 and L < 3, and therefore the con- 
tributing partial waves are, in a spectroscopic notation, 
3 P , 3 Pi, 3 P2~ 3 F 2 and X D 2 . 



Standard techniques [2l|, |43[ are now used to carry out 
the multipole expansion of the weak charge, /o(q), and 
current, j(q), operators. For muon capture on deuteron, 
we find 



(* L n n JJ °(p)\p(<l)\Msd)) 



(<f J »U,(q)|^(s d )) 



(^tn JJ '(.P)\jx(cO\^d(s d )) 



A>0 V2J + 1 

A>0 V2J + 1 

V2A + T i A (l5rf ' A - A| JJz) [-AAff gJ (g) + , 

A -^1 V + 1 



(2.13) 
(2.14) 
(2.15) 



r 



where A = ±1, and C^ SJ (q), 



E^ SJ (q) and 



M^ SJ (q) denote the reduced matrix elements (RMEs) 
of the Coulomb (C), longitudinal (L), transverse electric 
(E) and transverse magnetic (M) multipole operators, 
as defined in Ref. [2l|. Since the weak charge and cur- 



rent operators have scalar/polar- vector (V) and pseudo- 
scalar/axial- vector (A) components, each multipole con- 
sists of the sum of V and A terms, having opposite parity 
under space inversions |21| . The contributing multipoles 
for the S-, P-, and D-channels considered in the present 
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work are listed in Table U where the superscripts LSJ 
have been dropped. 

In the case of muon capture on 3 He, explicit expres- 
sions for the multipole operators can be found in Ref. [13] ■ 
Here we only note that parity and angular momentum 
selection rules restrict the contributing RMEs to Cb(V), 
d(A), L (V), L ± (A), E 1 {A), and M 1 (V) ; 

The total capture rate for the two reactions under con- 
sideration is then defined as 



dT = 2tt6{AE)\T w \ 2 x (phasespace) , 



(2.16) 



where 6(AE) is the energy-conserving (5-function, and 
the phase space is dp dk u /(2ir) 6 for reaction (jl.ll) and 
just <ik„/(27r) 3 for reaction (|1.2|) . The following notation 
has been introduced: (i) for muon capture on deuteron 



<-W\ 



BiB 2 hv f z 



(2.17) 



and the initial hyperfine state has been fixed to be / 
1/2; (ii) for muon capture on 3 He 



(2.18) 



',K ffz 



and the factor 1/4 follows from assigning the same prob- 
ability to all different hyperfine states. 



After carrying out the spin sums, the total rate for 
muon capture on 3 He reads (HI 



Gt Ei\\- 



E v 



|Kr \Co(V) - L„(V0r + \Cx(A) - Li(A)Y + \Mx{V) - E 1 {A)Y 



(2.19) 



with E v given by 



III. NUCLEAR WAVE FUNCTIONS 



E v = 



(m A 



'3H 



2(to m + mane) 



(2.20) 



In the case of muon capture on deuteron, the differential 
rate reads 



dT L 
dp 



where 



= E, 



1 - 



(m M + m d ) 



, p 2 dp 

8vr 4 



E v = 



(m M + m d ) 2 - -iml - Ap 2 
2(to m + m d ) 



(2.21) 



(2.22) 



In Eqs. (|2.19|) - (|2.22|) . m^, m n , m,d, ma^i, ™> 3 Tie are the 
muon, neutron, deuteron, 3 H and 3 He masses. The in- 
tegration over p in Eq. (|2.21l) is performed numerically 
using Gauss-Legendre points, and a limited number of 
them, of the order of 10, is necessary to achieve conver- 
gence to better than 1 part in 10 3 . In order to calculate 
the total capture rate T D , the differential capture rate 
is plotted versus p, and numerically integrated. Usually, 
about 30 points in p are enough for this integration in 
each partial wave. The differential cross section for each 
contributing partial wave will be shown in Sec. [V] 



Bound and continuum wave functions for both 
two- and three-nucleon systems are obtained with the 
hyperspherical-harmonics (HH) expansion method. This 
method, as implemented in the case of A = 3 systems, 
has been reviewed in considerable detail in a series of re- 
cent publications [H, Hf| EH • We will discuss it here in 
the context of A = 2 systems, for which, of course, wave 



Partial wave Contributing multipoles 



3 Pi 



3p 2 _3j 



l D 2 



Ci{A), Li(A), 
d(V), Li (V), 
C (A), 
Ci(V), Li(V), 
C a (A), L 2 (A), 
d(V), L X {V), 
C a (A), L 2 (A), 
C 3 (V), L 3 (V). 
Ch(A), Li(A), 
C 2 (V), L 2 (V), 
C a (A), L Z {A), 



Ei(A), 
Ei(V)- 
L (A), 
Ei(V), 

MA), 
Ei(V), 
E 2 (A), 
E 3 (V). 
Ei {A), 
E 2 (V), 
E 3 (A), 



Afi(V) 
Mi (A) 

M^A), 
M 2 (V) 
Mi(A), 
M 2 (V), 
M 3 (A) 
Mi(V), 
M 2 {A), 
M 3 (V) 



TABLE I: Contributing multipoles in muon capture on 
deuteron, for all the nn partial waves with J < 2 and L < 3. 
The spectroscopic notation is used. See text for further ex- 
planations. 
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functions could have been obtained by direct solution of 
the Schrodinger equation. 

The nuclear wave function for a bound system of total 
angular momentum JJ Z can be generally written as 



\* JJt ) = Y, c i>K J ") . 



(3.1) 



where \^j l Jz ) is a complete set of states, and u is an 
index denoting the set of quantum numbers necessary to 
completely specify the basis elements. The coefficients of 
the expansion can be calculated by using the Rayleigh- 
Ritz variational principle, which states that 



(5 c V JJz \ H-E\V 



, 



(3.2) 



where 8 c ^ JJz indicates the variation of ty JJ * for arbi- 
trary infinitesimal changes of the linear coefficients c M . 
The problem of determining c M and the energy E is then 
reduced to a generalized eigenvalue problem, 



"£(* J /*\H-E\* J /*)c^=0 



(3.3) 



In the case of the deuteron, we first define 

^LSJJ; (r) = F LS j{r)y LS .j.j z (r) , (3.4) 

where 

yLSJjAi) = \?L{i)®Xs]jj. , (3-5) 

and r is the relative position vector of the two nucleons. 
Obviously, for the deuteron, L — or 2, S = 1, J = 1, 
and T = (the isospin state r\TTz with T = T z = has 
been dropped for brevity). The radial functions Flsj(t) 
are then conveniently expanded on a basis of Laguerre 
polynomials as 



FhS.j(r) = 2J c LS j,i fi(r) , 



(3.6) 



with 



fi(r) 



(1 + 2) 



7 3 / 2 Lf\ ir ) e-^ 2 t (3.7) 



where the parameter 7 is variationally optimized (7 is in 
the range of 3-4 fin" 1 for the AV18 and N3LO poten- 
tials) . The complete wave function is then reconstructed 
as 



E E c ^v/iM;w*(*) 



L=0,2 I 



(3.8) 



and the subscript fj, denotes the set of quantum numbers 
[L,S= 1, J= 1,1}. 



The nn continuum wave function is written as 



^jLSJJ z ^jLSJJz 



LSJJz 



(3.9) 



where ^^ SJJz describes the system in the region where 
the two neutrons are close to each other and their mu- 
tual interactions are strong, while <J/^ SJJi describes their 
relative motion in the asymptotic region. The function 
, which vanishes in the limit of large separations, 



LSJJ, 

c 



is expanded on Laguerre polynomials as before for the 
case of the deuteron, see Eq. (|3.1|) . while the function 



^LSJJz j g ^ e appropriate asymptotic solution in chan- 
nel I.S.I. 



<T. LS J J z 

* A 



where 



with 



E 

US' 



(3.10) 



sfisjj =RL I ipr)y LSJJz {i), (3.11) 



■OLipr) 



Rl(pr) = p L+1 f R (r)n L ( P r) , 



(3.12) 



p being the magnitude of the relative momentum. The 
functions Ji(pr) and n^(pr) are the regular and ir- 
regular spherical Bessel functions, and ]r(t) = [1 — 
exp(— br)] 2L+1 has been introduced to regularize n^ipr) 
at small values of r. The trial parameter b is taken as 
6 = 0.25 fin" 1 . 

The matrix elements T^-ls l'S'(p) an( ^ the coefficients 

• • ' LSJJ 

entering the expansion of ^ c z are determined ap- 
plying the Kohn variational principle [13], stating that 
the functional 



k j ls , l , s ,(p) 



2u 



L'S'JJz 



H - 



2//. 



LSJJz 



, (3.13) 



is stationary with respect to variations of the trial pa- 
rameters in \$> LSJJ *. Here u is the reduced mass of the 
nn system. Performing the variation, a system of lin- 
ear inhomogeneous equations for c M and a set of alge- 
braic equations for TZ'l s L , s ,(p) are derived, and solved 

by standard techniques. From the It-is l<S'(p) ,s i phase 
shifts, mixing angles, and scattering lengths are easily 
obtained. 

We list in Tables [TT] and Mil binding energies and scat- 
tering lengths calculated with the Hamiltonian models 
considered in the present work. Note that A = 3 wave 
functions retain both T = 1/2 and T = 3/2 contribu- 
tions, T being the total isospin quantum number. The 
experimental data are, in general, quite well reproduced. 
It should be noted that, by using only two-nucleon in- 
teraction, the triton and 3 He binding energies are 7.624 
MeV and 6.925 MeV with the AV18, and 7.854 MeV and 
7.128 MeV with the N3LO, respectively. 
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AV18 N3LO 



Exp. 



B d (MeV) 2.22457 2.22456 2.224574(9) 
Pd (%) 5.76 4.51 

(fm) -18.487 -18.900 -18.9(4) 

np (fm) -23.732 -23.732 -23.740(20) 
5.412 5.417 5.419(7) 



1 



(fn 



TABLE II: Deuteron binding energy B d (in MeV) and D-state 
probability (in %), nn and singlet and triplet np scattering 
lengths (in fm), calculated with the two-nucleon potentials 
AV18 and N3LO. The experimental results are from Ref. [39^ . 



Finally, we define 



—LSJJ 



US' 



(3.14) 



LSJJ LSJJ 

so that \P ' z has unit flux. The function if? z enters 



in Eq. (|2.12[) and in the expression for differential muon 
capture rate on deuteron. 

AV18/UIX N3LO/N2LO Exp. 



B 3li (MeV) 


8.479 


8.474 


8.482 


B 3Hc (MeV) 


7.750 


7.733 


7.718 


2 a nd (fm) 


0.590 


0.675 


0.645(3)(7) 


4 a„d (fm) 


6.343 


6.342 


6.35(2) 



TABLE III: Triton and 3 He binding energies _B3 H and S3 Ho 
(in MeV), and nd doublet and quartet scattering lengths (in 
fm), calculated with the two- and three- nucleon potentials 
AV18/UIX and N3LO/N2LO. The experimental results are 
from Ref. [H. 



IV. WEAK TRANSITION OPERATOR 

The nuclear weak charge and current operators consist 
of polar- and axial-vector components. In the present 
work, we consider two different models, both of which 
have been used in studies of weak pp and hep capture re- 
actions in the energy regime relevant to astrophysics [20l — 
I22I ]. The first model has been developed within the 
so-called "Standard Nuclear Physics Approach" (SNPA) 
and has been applied also to study weak transitions for 
A = 6 and 7 nuclei [48| . and magnetic moments and 
Ml widths of nuclei with A < 7 [Igj, It will be dis- 
cussed in Sec. IIV Al In the second model, the nuclear 
weak transition operators have been derived in heavy- 
baryon chiral perturbation theory (HBChPT), carrying 
out the expansion up to next-to-next-to-next-to leading 
order (N3LO) [22l l5l|. This same model has been up- 
dated and most recently used in the electromagnetic sec- 
tor to study the three-nucleon magnetic moments and 



np and nd radiative capture reactions in Ref. 
review it in Sec. IIVBI 



We 



A. The "Standard Nuclear Physics Approach" 

The one-body axial charge and current operators have 
the standard expressions [21fl obtained from the nonrela- 
tivistic reduction of the covariant single-nucleon current, 
and include terms proportional to 1/m 2 , m being the nu- 
cleon mass. The induced pseudoscalar contributions are 
retained both in the axial current and charge operators, 
and are given by 



2 m m h 

9Ps{ql) 
2 m nil, 



(4.1) 

T i ,_q(ir i .q)e i * r4 , (4.2) 



where 



(n, 



IT, 



'■y 



)/2 



(4.3) 



and q is the three-momentum transfer, m M the muon 
mass, E v the neutrino energy, g 2 . the squared four- 
momentum transfer. Note that, since in the present case 



q is not negligible (q 



rn 



t ), axial and induced pseu- 



doscalar form factors need be included, 
of Ref. [34J , they are taken as 



In the notation 



9A{ql) 

gps{ql) 



9A 



(l + ??/Ai) 2 ' 

■9A(q a ) 



For 



(4.4) 
(4.5) 

the axial- vector coupling constant 17^4, two val- 
ues have been adopted in this study: the first one, 
#,4=1. 2654(42), is taken from Ref. [53, and has been used 
widely in studies of weak processes [2ll. I22I l3~il |48| . The 
second one, <M=1- 2695(29), is the latest determination 
quoted by the Particle Data Group (PDG) [Ej|- The two 
values for qa are consistent with each other. In order to 
compare with the results of the aforementioned studies, 
we have adopted the earlier determination for qa- On the 
other hand, to estimate the theoretical uncertainty aris- 
ing from this source, the calculation has been carried out 
also using (^=1.2695(29) in one specific case. As it will 
be shown below, the central value and the theoretical un- 
certainties of the considered observables are comparable 
in both cases. 

The value for the cutoff mass Aa used in this work is 
1 GeV/c 2 , as in Ref. [34|. It is obtained from an analysis 
of pion electroproduction data 1541 a nd measurements of 
the reaction + p — > n + u + 55]. Since here q%. << 
A A , an uncertainty of few % on A^ is expected to affect 
gA(q&) at the percent level. The ^-dependence of gps 
is obtained in accordance with the partially-conserved- 
axial-current (PCAC) hypothesis, by assuming pion-pole 
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dominance and the Goldberger-Treiman relation [43j, |56| . 
In Eq. (|4.5I) . is the pion mass. 

The two-body weak axial-charge operator includes 
a pion-range term, which follows from soft-pion theo- 
rem and current algebra arguments [H?], HH , and short- 
range terms, associated with scalar- and vector- meson 
exchanges. The latter are obtained consistently with 
the two-nucleon interaction model, following a proce- 
dure [59[ similar to that used to derive the corresponding 
weak vector-current operators [2l|. The two-body axial 
charge operator due to 7V-to-A-transition has also been 
included [2l], [34| , although its contribution is found to be 
very small. 

The two-body axial current operators can be divided 
in two classes: the operators of the first class are derived 
from 7T- and p-meson exchanges and the /?7r-transition 
mechanism. These mesonic operators, first obtained in a 
systematic way in Ref. [60] , give rather small contribu- 
tions [21] . The operators in the second class are those 
that give the dominant two- body contributions, and are 
due to A-isobar excitation [2(| [2l[ . We review them here 
briefly. The A-to-A-transition axial current is written as 
(in the notation of Ref. [2lj ]) 



j z (1) (q;iV^ A, A) 



3a(<7 ct )S 1 + .J* q(S, • q) 



2 777 777. 



x e 



(4.6) 



where Si and Tj are spin- and isospin-transition opera- 
tors, which convert a nucleon into a A-isobar. The in- 
duced pseudoscalar contribution has been obtained from 
a non-relativistic reduction of the covariant iV-to-A- 
transition axial current [Hi]. 

The axial and pseudoscalar form factors g* A and g* PS 
are parameterized as 



gpsiil) 



= RA9A{q a ) 

2 77l„ 777 



gM) 



(4.7) 



with <7a(<2£) given in Eq. (|4.4jl . The parameter Ra 
is adjusted to reproduce the experimental value of 
the Gamow- Teller matrix element in tritium /3-decay 
(GT EXP ), while the g 2 -dependence of g* PS is again ob- 
tained by assuming pion-pole dominance and PC AC l43l . 
[HI]. The value for GT EXP , estimated in Ref. [2(3], 
was 0.957(3). This value was determined by assuming 
&4=1.2654(42), (F) 2 =0.9987, where (F) is the reduced 
matrix element of the Fermi operator yv t ,;._, and the 
triton half-life /T 1/2 is (1134.6±3.1) s [6l|].' We adopt 
it also in the present work, except for the case in which 
g A is taken from the PDG, 5^=1-2695(29). We then ex- 
tract GT EXP = 0.955(2), corresponding to (F) 2 =0.99926 
and /T 1/2 =(1132.1±4.3) s. The new value of (F) 2 dif- 
fers by less than 0.1 % from the older, presumably due 
to the higher accuracy of the present trinucleon wave 
functions. The new value of /I\/ 2 has been obtained 
by averaging the previous value of fTi/ 2 with the new 



one of Ref. (62[, (1129. 6±3) s, and summing the errors in 
quadrature. The values for Ra determined in the present 
study using trinucleon wave functions corresponding to 
the AV18/UIX are ^=1-21(9), when GT EXP = 0.957(3) 
and 0a=1<2654(42), and i? A =1.13(6), when GT EXP = 
0.955(2) and (^=1-2695(29). The experimental error on 
GT EXP is responsible for the 5-8 % uncertainty in Ra- 

It is important to note that the value of Ra de- 
pends on how the A-isobar degrees of freedom are 
treated. In the present work, as in Ref. 34], the two- 
body A-excitation axial operator is derived in the static 
A approximation, using first-order perturbation theory 
(PT). This approach is considerably simpler than that 
adopted in Ref. [2l[, where the A degrees of freedom 
were treated non-perturbatively, within the so-called 
transition-correlation operator (TCO) approach, by re- 
taining them explicitly in the nuclear wave functions [63j . 
The results for Ra obtained within the two schemes dif- 
fer by more than a factor of 2 [2l|. However, the re- 
sults for the observables calculated consistently within 
the two different approaches are typically within 1 % of 
each other. Finally, because of the procedure adopted 
to determine Ra, the coupling constant g* A of Eq. (14. 7f> 
cannot be naively interpreted as the AA axial coupling 
constant. The excitation of additional resonances and 
their associated contributions will contaminate the value 

oig* A - 

The weak vector charge and current operators are con- 
structed from the isovector part of the electromagnetic 
current, in accordance with the conserved- vector-current 
(CVC) hypothesis. The weak charge operator includes 
the non-relativistic one-body term and the relativistic 
spin-orbit and Darwin- Foldy contributions, and is ob- 
tained from the corresponding isovector electromagnetic 
operator, as listed in Ref. [64j, by replacing 



n, z /2 



(4.8) 



Two-body contributions, arising from n- and p- meson ex- 
change mechanisms, are also included. Their expressions 
are listed in Ref. (64j, with the substitution 



and 



(Ti X Tj), -> (Ti X Tj )- , 



(Tj X Tj)- = (Tj X Tj) x - i(Ti X Tj) y 



(4.9) 



(4.10) 



Electromagnetic form factors are also included, and avail- 
able parametrizations for them all provide excellent fits 
of the experimental data at the low momentum transfer 
of interest here. 

The weak vector current operator retains the one- 
body operator, two-body "model-independent" (MI) and 
"model-dependent" (MD) terms, and three-body terms. 
The MI two-body currents are obtained from the two- 
nucleon interaction, and by construction satisfy current 
conservation with it. In the present work, we include 
the leading two-body "7r-like" and "p-like" operators, ob- 
tained from the isospin-dependent central, spin-spin and 
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tensor nucleon-nucleon interactions. On the other hand, 
we have neglected the additional two-body currents aris- 
ing from non-static (momentum-dependent) interactions, 
since these currents are short-ranged, and numerically far 
less important than those considered here [49j |. The MD 
currents are purely transverse, and therefore cannot be 
directly linked to the underlying two-nucleon interaction. 
The present calculation includes the isovector currents 
associated with the ujttj transition mechanism [64]], and 
the excitation of A isobars. The former contributions are 
numerically negligible, while the latter are important in 
order to reproduce the three-nucleon magnetic moments 
and elastic form factors [4{|. The contributions of the 
(MD) A-isobar currents have been calculated in this case 
with the TCO method [&3l, and explicit expressions for 
these are listed in Ref. [65|. Again, the substitution of 
Eq. (|4~T0t is used. 

A three-nucleon interaction requires a corresponding 
three-body current. The latter was first derived in 
Ref. [49j from a three-nucleon interaction consisting of 
a dominant two-pion-exchange component, such as the 
UIX model adopted in the present work. The charge- 
changing three-body (weak vector) current is obtained 
applying CVC to the operators listed in Ref. [49j], i.e. 



■ Tb, 



j§L(q; 7) 



Ea*&4 3 iL(q; V) 



(4.11) 

(3) 

where j-^ 2 (q;7) are the isovector three-body electro- 
magnetic currents, and a,b — x, y, z are isospin Cartesian 
components. 

We conclude by emphasizing that the model for the 
weak transition operator is the same as that of Ref. [13] , 
but for two differences relative to its vector component: 
(i) the MD two-body currents due to A-isobar degrees 
of freedom are treated non perturbatively with the TCO 
method, rather than in first order PT, and (ii) three- 
body terms are also included - they were neglected in 
Ref. [2H . As shown in Table IIV1 the present model for 
the electromagnetic current provides an excellent descrip- 
tion of the trinucleon magnetic moments, in particular 
of their isovector contribution, to better than 1 % (row 
labeled with "FULL"). This gives us confidence in the 
accuracy of the corresponding weak vector currents. We 
also note that: (i) three-body contributions are impor- 
tant to achieve this level of agreement between theory 
and experiment; (ii) the results obtained with the model 
of Ref. [13, labeled by "MD/A-PT", are at variance with 
data at the 3 % level; (iii) the difference between the 
present "FULL" results and those of Ref. [5i| can be 
traced back to the fact that here we have used slightly 
more accurate wave functions, with a careful treatment of 
T = 3/2 components. However, we should note that, in 
contrast to Ref. GPL MI two-body contributions arising 



from non-static components of the two-nucleon interac- 
tion have been ignored. They amount to a negligible 0.1 
% correction. 



3 He 



IA 2.5745 -1.7634 
MI 2.8954 -2.0790 
MD/A-PT 3.0260 -2.2068 
MD/A-TCO 2.9337 -2.1079 
FULL 2.9525 -2.1299 
Ref. [49] 2.953 -2.125 



Exp. 



2.9790 -2.1276 



TABLE IV: Triton and 3 He magnetic moments, in nu- 
clear magnetons, calculated with the AV18/UIX Hamiltonian 
model and compared with the experimental data. The re- 
sults labeled "I A" are obtained with single-nucleon current, 
while those labeled "MI" are obtained by including in ad- 
dition the model-independent two-body contributions. The 
results labeled "MD/A-PT" and "MD/A-TCO" include also 
the model-dependent contributions, with the A-isobar degrees 
of freedom treated in perturbation theory or within the TCO 
approach, respectively. Lastly, the results labeled "FULL" 
retain three-body contributions. Also shown are the results 
of Ref. [H. 



B. The chiral effective field theory approach 

ThevEFT weak transition operator is taken from 
Rcfs. [22j and 142], where it was derived in covariant 
perturbation theory based on the heavy-baryon formu- 
lation of chiral Lagrangians by retaining corrections up 
to N3LO. It was recently used in Ref. [66| to study Ml 
electromagnetic transitions. We review its main features 
here. 

The vector and axial- vector one-body operators are the 
same as those obtained within the SNPA, and described 
in Sec. IIV Al These one-body operators are also listed 
in Eq. (17) of Ref. [2^], except that we also include cor- 
rections in the vector current, which arise when the non- 
relativistic reduction is carried out to next-to-leading or- 
der (proportional to 1/m 3 ). The resulting operator is 
given (in momentum space) by 

Ji '(q;V;=Ji (q;^)+j, (q;V) , (4.12) 

where j^^^q; V) is the standard leading order term 



Ji 



(1)NB. 



(q;V) 



GMjk + iGM^xq 



V / 2n 



and j- 1)iiC (q;U) is 



(4.13) 
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.(1)RC 



4m 3 



G v E {ql) k 2 (2k + iff, x q) + i[G&(<£) - G£(<£)](2k • q<7, X k + k X qtr* • k) 



In the above expressions, k = (p- +Pi)/2, q = p- p^ p- 
and Pi being the momenta of the outgoing and ingoing 
nucleons, respectively, and G^(ql) and G\ lI {q^) denote 
the isovector combinations of the nucleon electric and 
magnetic form factors [64j. Note that Eqs. (14.121) (14. 14|) 
can also be obtained from Eq. (11) of Ref. [42], with the 
substitution (14.81) required by CVC. 

The axial two-body charge and currents are from 
Ref. [22j. In particular, the axial charge operator is 
that derived originally in Ref. [13] ■ In the SNPA, ad- 
ditional contributions arc considered, which, in a %EFT 
context, are expected to appear at higher orders. The 
two-body axial current operator consists of two contribu- 
tions: a one-pion exchange term and a (non-derivative) 
two-nucleon contact-term. The explicit expressions for 
these terms can be found in Ref. (22[. While the cou- 
pling constants which appear in the one-pion exchange 
term are fixed by ttN data, the low-energy constant da, 
determining the strength of the contact-term, has been 
fixed by reproducing GT EXP . The value ^=1.2654(42) 
has been used. The values for da are presented in Ta- 
ble [V] and discussed below. 

The two-body vector currents are decomposed into 
four terms (42J: the soft-one-pion-exchange (l7r) term, 
vertex corrections to the one-pion exchange (l7rC), the 
two-pion exchange (27r) , and a contact-term contribution. 
Their explicit expressions can be found in Eqs. (12)-(18) 
of Ref. [42] , with the replacements (I4.3[) and (|4.10p in the 
isospin operators. All the Itt, InC and 27r contributions 
contain low-energy constants estimated in Ref. [671 ] . us- 
ing resonance saturation arguments. The contact-term 
electromagnetic contribution is given by 

4 2)CT (q;7) = -=j-e 1 * R qx[ 5 4 S ( < r i + < r i ) 
J Ira 

+ 94v{ti x Tj) z <Ti x tTj)8\{r) . (4.15) 

The function S\(r), as well as the Yukawa functions 
which appear in the In, InC and 2-7T operators and in the 
two-body axial current terms, are obtained by performing 
the Fourier transform from momentum- to coordinate- 
space with a Gaussian regulator characterized by a cut- 
off A. This cutoff determines the momentum scale below 
which these EFT currents are expected to be valid, i.e., 
A=500-800 MeV [H. The explicit expression of S A {r) 
in Eq. ([4~T5|) is 

*A(r) = /(^3 e " fc2/A2eikr - (4- 16 ) 

The coefficients g±s and g±v of Eq. (14. 15)) are fixed 
to reproduce the experimental values of triton and 3 He 



(4.14) 

I 

magnetic moments, for each nuclear interaction and cut- 
off value. This procedure is similar to that used to fix 
the strength da of the contact-term in the two-body ax- 
ial current discussed above. The results are presented in 
Table [V] A few comments are in order: (i) while values 
of both <?4s and g^y are presented, only g$y is relevant 
in the present work, since CVC relates the isovector elec- 
tromagnetic current to the weak vector current, (ii) The 
uncertainty on g^s and g^v is not due to the experimental 
errors on the triton and 3 He magnetic moments, which 
are in fact negligible, but rather to numerics. We have 
used a random walk consisting of 1.6 million configura- 
tions, in order to reduce the numerical uncertainty on 
A/i, the difference between the experimental magnetic 
moments and the result obtained without the contact 
contributions, to less than 1 %. In contrast, the exper- 
imental error on GT EXP is primarily responsible for the 
uncertainty in du. (iii) The values reported for du are dif- 
ferent, but consistent within the error, with those listed 
in Table II of Ref. [22j . This is due to differences be- 
tween t he p resent 3 H and 3 He wave functions and those 
of Ref. [22|] - we have already commented on this in the 
previous section, (iv) The g^s and t/4y values are rather 
different from those listed in Ref. [66| . We observe that 
<74S and g 4 y are fixed by fitting very small quantities. In 
the AV18/UIX case with A = 600 MeV, for instance, the 
value A/i is 0.0461(4) for triton, and -0.0211(4) for 3 He, 
where the uncertainties are statistical errors due to the 
Monte Carlo integrations. Consequently, the resulting 
values will be sensitive to several factors, including nu- 
merics. However, it is worthwhile pointing out that the 
isovector contact contribution to the muon capture rates 
under consideration turns out to be negligible. Finally, 
we should point out that the present %EFT model for the 
weak vector current operator differs in some of its two- 
pion exchange parts from that obtained in time-ordered 
perturbation theory by some of the present authors in 
Ref. |68[. The origins of these differences have been dis- 
cussed in Refs. 68, 69]. However, for consistency with 
the calculations of the pp and hep reactions of Ref. [22j , 
we have chosen to use the x^FT model illustrated above. 
We do not expect these differences to be numerically sig- 
nificant for the processes under consideration. 



V. RESULTS 

The results for the total rates of muon capture on 
deuteron and 3 He are presented in the following two sub- 
sections. 



11 



A (MeV) gis gw d R 





500 


0.69(1) 2.065(6) 


0.97(7) 


AV18/UIX 


600 


0.55(1) 0.793(6) 


1.75(8) 




800 


0.25(2) -1.07(1) 


3.89(10) 


N3LO/N2LO 


600 


0.11(1) 3.124(6) 


1.00(9) 



TABLE V: The LECs g^s and g&v associated with the 
isoscalar and isovector contact terms in the electromagnetic 
current (see Eq. (|4.15|l 'l. and the LEC cLr, of the two-body 
axial-current contact term, calculated for three values of the 
cutoff A with triton and 3 He wave functions obtained from the 
AV18/UIX model. For A = 600 MeV, also the N3LO/N2LO 
model is also used. 



SNPA (AV18) ]So 3 P 3 Pi 3 P 2 1 D 2 3 F 2 Total 

3A=1.2654(42) 246.6(7) 20.1 46.7 71.6 4.5 0.9 390.4(7) 

gA=1.2695(29) 246.8(5) 20.1 46.8 71.8 4.5 0.9 390.9(7) 

EFT* (AV18) ]Sq 3 P 3 Pi 3 P 2 1 D 2 S F 2 Total 

A = 500 MeV 250.0(8) 19.9 46.2 71.2 4.5 0.9 392.7(8) 

A = 600 MeV 250.0(8) 19.8 46.3 71.1 4.5 0.9 392.6(8) 

A = 800 MeV 249.7(7) 19.8 46.4 71.1 4.5 0.9 392.4(7) 

EFT* (N3LO) ]So 3 Po 3 Pi 3 P 2 1 D 2 Z F 2 Total 

A = 600 MeV 250.5(7) 19.9 46.4 71.5 4.4 0.9 393.6(7) 

TABLE VI: Total rate for muon capture on deuteron, in the 
doublet initial hyperfine state, in s _1 . The different partial 
wave contributions are indicated. The numbers among paren- 
theses indicate the theoretical uncertainty arising from the 
adopted fitting procedures, as explained in Sec. IIVI Such un- 
certainty is not indicated when less than 0.1 s _1 . The label 
"SNPA" indicates that the results have been obtained using 
the model for the weak transition operator of Sec. IIV Al while 
the label "EFT*" indicates that the model of Sec. lIVBl is used. 
The AV18 and N3LO interactions have been used to calculate 
the deuteron and nn wave functions. 

A. Muon capture on deuteron 

In a partial wave expansion of the final nn state, all 
channels with total angular momentum J < 2 and rel- 
ative orbital angular momentum L < 3 have been in- 
cluded, i.e., 1 S , 3 P , 3 Pi, 3 P2, 1 D 2 and 3 F 2 . Partial 
waves of higher order contribute less ~ 0.5 % to the rate. 
Indeed, the 3 F 2 contribution turns out to be already be- 
low this level. 



We present in Table IVII the results for the total cap- 
ture rate in the doublet hyperfine state. Both models 
for the nuclear weak transition operator presented in 
Sec. HVAl and Sec. HVBl have been used, labeled SNPA 
and EFT*, respectively. The nuclear wave functions have 
been calculated with the AV18 [H or the N3LO (H two- 
nucleon interaction. The label EFT* is used to denote 
the results of calculations in which the matrix elements of 
xEFT weak operators are evaluated between wave func- 
tions corresponding to both conventional and chiral po- 



tentials. The first approach (based on a conventional 
potential) is often referred to in the literature as the "hy- 
brid" approach. The second, using chiral potentials and 
currents, is in principle a full-fledged %EFT calculation, 
except that these potentials and currents have not (yet) 
been derived consistently at the same order in the low- 
momentum scale. For this reason, we characterize the 
corresponding results with the EFT* label. 

Within each Hamiltonian model, the parameters 
present in the SNPA and EFT* axial current models 
have been fitted to reproduce GT EXP in tritium /3-decay, 
as discussed in Sec. IIVI Furthermore, the LECs in the 
EFT* weak vector current have been fitted to reproduce 
the A = 3 magnetic moments. The three-nucleon wave 
functions have been generated, in this fitting procedure, 
from two- and three-nucleon interactions, either AV18 
and UIX [H, or N3LO and N2LO 40]. Inspection of 
the table shows that the 1 Sq contribution is the leading 
one, but L > 1 contributions are significant and account 
for ~ 37 % of the total rate. By comparison between 
the first and second row of the table, we conclude that 
there is no difference in the results, within uncertainties, 
when the older value for gA, <M=L2654(42), or the most 
recent one, (^=1.2695(29), is used. This reflects the fact 
that the factor Ra, see Sec. IIV A\ has been constrained 
by GT EXP in both cases. The model dependence due to 
interactions, currents, and the cutoff A (present in the 
%EFT version of these currents) is at the 1 % level in 
the total rate, and hence very weak. It is a bit larger, 2 
%, in the 1 Sg channel, for which two-body current con- 
tributions are larger, see Tables IVlIl and [Villi below. In 
conclusion, a total capture rate in the range 

T D = (389.7 - 394.3) s" 1 , (5.1) 

can be conservatively ascribed to reaction (jl.ip . This re- 
sult is in agreement with the measurements of Refs. [1, H 
,8j , but not with that of Ref. Mj. The differences with the 
theoretical results of Refs. [R^ 4l8| _are also very small. In 
particular, the authors of RefipJ] find r D ( 1 S' ) = 245 
s — 1 , in contrast to T D ( 1 So) = 250 s _1 reported here. We 
have explicitly verified that this difference is mainly due 
to the inclusion in the present work of mesonic two-body 
contributions in the weak vector current beyond the soft 
one-pion exchange term discussed in Sec. IIVBI On the 
other hand, the results of Refs. [HI, [l{| are significantly 
larger than those listed here, presumably because these 
authors have not constrained their weak current to repro- 
duce GT EXP and the isovector magnetic moment of the 
trinucleons. We observe that our approach also provides 
a value for muon capture on 3 He in excellent agreement 
with the experimental datum. 

For future reference, we show in Fig. Q] the differential 
capture rate dT D /dp, defined in Eq. (|2.21|) . as function of 
the nn relative momentum p, calculated with the SNPA 
(AV18) model. Integrating each partial wave contribu- 
tion leads to the values listed in the first row of Table IVT1 




The cumulative contributions to the leading 1 So and 3 P 2 
capture rates of various terms in the nuclear current oper- 
ator are given in Table lVIII for the SNPA and lVIIII for the 
EFT*. Table ETI] shows that: (i) the (one-body) induced 
pseudoscalar contribution is significant, as expected, and 
reduces the rate of - 20 % (~ I 3 %) in x 5 ( 3 P 2 ) cap- 
ture; (ii) the mesonic contributions are very small; (iii) 
the A-isobar contributions are significant only in S-w&ve 
capture, and differences in the vector two-body current 
originating from treating the A-isobar degrees of free- 
dom either perturbatively or within the TCO approach 
are negligible; (iv) the A-isobar induced pseudoscalar ax- 
ial current contribution of Eqs. (|4.6p and (|4.7[) is of the 
order of I s™ 1 for 1 So capture. 

From inspection of Table [VTlTT we can conclude that: 
(i) the one-body contributions (labeled "IA+RC" and 
"PS") are slightly different from those reported in the 
SNPA. While the EFT* one-body axial charge and cur- 
rent, and vector charge, operators are the same as in 
the SNPA, the EFT* one-body vector current includes 



the relativistic corrections of Eq. (14. 14[) . ignored in the 
SNPA. On the other hand, the differences between SNPA 
and EFT* one-body contributions remain well below the 
1 % level, and therefore within the theoretical uncer- 
tainty. However, it would be interesting to study the 
contributions of these relativistic corrections to the many 
electromagnetic observables analyzed within the SNPA 
over the years [49|, 15011 - A first step in this direction has 
been done in Ref. [70J. (ii) The A dependence is very 
weak, well below the 1 % level, (iii) The 3 P 2 term does 
not show any A dependence nor any sensitivity to the da 
and g4v fitting procedure. This is easily understood by 
observing that the contact terms in the vector and ax- 
ial currents are proportional to (Tj x Tj)^((Ti x crj) (see 
Ref. [22| and Eq. (14451) ). whose matrix elements vanish 
when calculated between the deuteron state and 3 Pj nn 
states. 

By comparison of the results listed in Tables IVIII 
and I VIIH we observe that the main source for the 2 % 
theoretical uncertainty in the Sq channel quoted above 
is the A-isobar induced pseudoscalar contribution: were 
it to be neglected, the SNPA and EFT* 1 S results would 
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IA+RC 293.0 82.6 

PS 233.2 71.6 

mesonic 237.5 71.4 

A-PT(A-w/o PS)+A-PT(V) 248.3(8) 71.6 
A-PT(A-w/o PS)+A-TCO(V) 247.7(8) 71.6 
A-PT(A-w PS)+A-TCO(V) 246.6(7) 71.6 

TABLE VII: Cumulative contributions to the total rate for 
muon capture on deuteron in s" 1 , when the model of Sec. lIV Al 
is used for the weak transition operator. The deuteron and 
nn wave functions are calculated using the AV18 interaction. 
Only the 1 So and 3 P 2 partial waves are considered, since they 
give the leading contributions. The labels "IA+RC", "PS", 
and "mesonic" indicate the results obtained by retaining the 
impulse approximation plus relativistic corrections, the in- 
duced pseudoscalar, and the purely mesonic contributions. 
The label "A-PT(yl-w/o PS)" is used to indicate that the 
axial charge and current contributions due to A-isobars are 
retained perturbatively, but the induced pseudoscalar con- 
tribution in the one-body A-current is not included, while 
the label "A-PT(^-w PS)" is used when the iV-to-A induced 
pseudoscalar contribution is also included. Finally, the la- 
bels 'A-PT(V)" and 'A-TCO(V)" are used when the vector 
current contributions due to A-isobars are treated perturba- 
tively and within the TCO scheme, respectively. The value 
(?a = 1-2654(42) is used. The numbers in parentheses indi- 
cate the theoretical errors arising from the fitting procedure 
adopted for Ra- This uncertainty is not shown if less than 
0.1 s" 1 . 

agree at the 1 % level. We observe that the strength of 
the associated operator, which is only included in the 
SNPA calculation and plays no role in tritium /3-decay 
(being proportional to the momentum transfer q), fol- 
lows from PCAC and pion-dominance arguments, see 

Sec. EES 

B. Muon capture on 3 He 

We present in Table lLXl thc results for the total capture 
rate Tq defined in Eq. (I2.19[) . The theoretical uncertain- 
ties due to the fitting procedure of Ra (SNPA) or cIr 
and gw (EFT*) are in parentheses. The results in the 
first two rows have been obtained in the SNPA, using the 
AV18/UIX Hamiltonian and the two available values for 
cja- These results are the same, since a change in qa is 
compensated by a corresponding one in Ra- The results 
of the next three rows have been obtained within the 
hybrid approach EFT* (AV18/UIX). They show a very 
weak A-dependence, and are in excellent agreement with 
those reported in SNPA. The result in the last row is with 
the EFT* (N3LO/N2LO). The EFT* (AV18/UIX) and 
EFT* (N3LO/N2LO) differ by 8 s" 1 , or less than 1 %. 
In view of this, we quote conservatively a total capture 
rate for reaction (11.21) in the range 

r = (1471 - 1497) s" 1 , (5.2) 



: 5o 3 P 2 







IA+RC 


292.5 82.5 


A = 


500 MeV 


PS 


232.7 71.6 






mesonic 


250.0(8) 71.2 






IA+RC 


292.5 82.5 


A = 


600 MeV 


PS 


232.7 71.6 






mesonic 


250.0(8) 71.1 






IA+RC 


292.5 82.5 


A = 


800 MeV 


PS 


232.7 71.6 






mesonic 


249.7(7) 71.1 



TABLE VIII: Same as T able IVIII but with the nuclear cur- 
rent model of Sec. ITVBl The labels "IA+RC", "PS", and 
"mesonic" indicate the results obtained by retaining the im- 
pulse approximation plus relativistic corrections, the induced 
pseudoscalar, and the mesonic contributions. The relativistic 
corrections to the one-body vector current operator include 
the additional terms of Eq. (|4.14[) . The numbers in paren- 
theses indicate the theoretical errors arising from the fitting 
procedure adopted for da and qav ■ This uncertainty is not 
shown if less than 1 s _1 . 



by keeping the lowest and upper bounds in the values of 
Table HI 



SNPA (AV18/UIX) 


To 


<M=1.2654(42) 


1486(8) 


5a=1.2695(29) 


1486(5) 


EFT* (AV18/UIX) 


Po 


A = 500 MeV 


1487(8) 


A = 600 MeV 


1488(9) 


A = 800 MeV 


1488(8) 


EFT* (N3LO/N2LO) 


To 


A = 600 MeV 


1480(9) 


EXP. 


1496(4) 



TABLE IX: Total rate for muon capture on 3 He, in s _1 . The 
numbers in parentheses indicate the theoretical uncertainties 
due to the adopted fitting procedure, see Sec. IIVI The la- 
bel "SNPA" ("EFT*") indicates that the results have been 
obtained using the model for the weak transition operator of 
Sec. HVAl fSec. lIVBX The triton and 3 He wave functions are 
obtained from the AV18/UIX and N3LO/N2LO Hamiltoni- 
ans. 



The contributions of the different components of the 
weak current and charge operators to the total rate and 
to the reduced matrix elements (RMEs) of the contribut- 
ing multipoles are reported in Table |X] for SNPA and IXII 
for EFT*. The HH wave functions have been calculated 
using the AV18/UIX Hamiltonian model. In TableE] R A 
has been fixed to its central value 1.21, and in Table IXfl 
we use A = 600 MeV and consequently — 1.75 and 
= 0.793 (see Table fV|) . The notation in these tables 
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is the same as that of Tables EED and IVTIll of Sec. IVAl 
with the only exception of the label "A-TCO(V)" , which 
here indicates that the two-body A-currents are treated 
with the TCO method, and that the three-body current, 
constructed consistently with the UIX three-nucleon in- 
teraction, are also included. In Table [Xj we observe 
that the induced pseudoscalar term gives a significant 
contribution to the L\(A) and C\(A) RMEs, although 
the latter is much smaller than L\(A) in magnitude. 
The mesonic contributions to C (V), Ci(A), L 1 (A), and 
Ei(A) are small, while they provide a 15 % correction 
to Mi(V), as expected (in the isovector magnetic mo- 
ments of the trinucleons, these mesonic currents give a 
15 % contribution relative to the one-body). Contribu- 
tions due to A-isobar excitations in the weak vector and 
axial currents are at the few % level. Finally, we note 
that the result for Tq is in excellent agreement with that 
of the earlier study of Ref. [34]. This agreement comes 
about because of two compensating effects: on the one 
hand, the Li(A) and E\(A) RMEs are, in magnitude, 
slightly larger, as consequence of the fact that the pa- 
rameter Ra is slightly larger here than in Ref. 34 1, be- 



and A (nucleon only) degrees of freedom. Indeed, the A- 
TCO(y) Mi(V) RME should be compared with M X {V) 
corresponding to A-PT(V) in Table Ixl 



cause of the more accurate wave functions employed. On 
the other hand, the M\(V) RME corresponding to A- 
TCO(V) is slightly smaller than calculated in Ref. [Hj], 
Mi(V) = 0.1355 in Ref. Q versus M X {V) = 0.1346. 
This last value results from dividing the value listed in 
Table IXl by the factor 1.022, arising from the normaliza- 
tion correction of the trinucleon wave functions due to 
the presence of explicit A-isobar degrees of freedom, i.e., 

N A = V(*A r +A|*A r +A)/(*Aronly|*7Vo„ly), where V N+A 

(^jVoniy) is the nuclear wave function with both nucleon 



In Table IXfl we observe that: (i) the one-body vector 
current contribution to the M\(V) RME is different than 
in SNPA. This is due to the presence of the additional 
rclativistic corrections of Eq. (|4.14|) . (ii) The mesonic 
contributions are significant for the Li(A), Ei(A) and 
Mi(V) RMEs, and bring their values closer to those in 
the SNPA. (hi) The lirC, 2tt, and contact terms in the 
mesonic vector current are important. If they were to 
be neglected, the total Mi(V) RME would be equal to 
0.1201 and consequently To = 1453 s _1 . (iv) The results 
of Table IXIl should be compared with those of Table 2 
of Ref. [33. We find significant differences in all the 
RMEs, both for the one-body contribution (here labeled 
"PS" and in Ref. [13 "IA") and the complete calculation. 
Only the results for the Mi (V) RME appear to be similar 
to each other, although it is unclear whether in Ref. [37j 
the vector l7rC, 2-7T and contact-term contributions are 
included. 

Finally, we observe that when the N3LO/N2LO Hamil- 
tonian model is used, C (V), C\(A), Li(A), Ei(A) and 
Mi(V) are -0.3288, 0.4130xl0" 2 , -0.2773, -0.5810 and 
0.1329, respectively, when the mesonic contributions are 
included. Comparing these results with the AV18/UIX 
ones of Table IXIl we see that all the RMEs are compa- 
rable at the 1 % level. This fact is reflected in the ~ 1 
% difference between the AV18/UIX and N3LO/N2LO 
results for Tq. Therefore, it does not seem possible to 
identify a particular source for this 1 % difference. 



r Cp(V) CijA) Li{A) ExjA) Mi(V) 

IA+RC 1530 -0.3287 0.7440xl0~ 2 -0.4056 -0.5516 0.1127 

PS 1316 0.3986xl0~ 2 -0.2589 

mesonic 1385 -0.3283 0.4024xl0~ 2 -0.2619 -0.5562 0.1315 

A-PT(t4-w/o PS)+A-PT(V / ) 1501 0.4287xl0" 2 -0.2811 -0.5821 0.1370 

A4 3 T(A-w/o PS)+A-TCO(V) 1493 0.1376 

A-PT(A-w PS)+A-TCO(V) 1486 -0.2742 

TABLE X: Cumulative contributions to the total rate To for muon capture on He, in s _1 , and to the reduced matrix elements 
(RMEs) Co(V), Ci(A), Li(A), E 1 (A) and M\(V) (see Sec. [TTJ), when the model of Sec. IIVAI is used for the weak transition 
operator (Ra = 1.21). The HH wave functions have been calculated using the AV18/UIX Hamiltonian. Notation is the same 
as in Table [VTT1 except for "A-TCO(V r )" I which here is used to indicate that the A-isobar contributions are treated with the 
TCO method, and, in addition, three-body weak vector current contributions are included. Note that Co(V) is purely real, 
while the other RMEs are purely imaginary. 



VI. SUMMARY AND CONCLUSIONS 



Total rates for muon capture on deuteron and 3 He 
have been calculated within a consistent approach, based 



15 



r Cq(V) Ci(A) L!(A) Ei{A) Mi(V) 

IA+RC 1517 -0.3287 0.7440xl0 -2 -0.4056 -0.5516 0.1082 

PS 1303 0.3986xl0~ 2 -0.2589 

mesonic 1488 0.3978xl0~ 2 -0.2810 -0.5833 0.1317 

TABLE XI: Same as Table [X] but with the nuclear current 
model of Sec. IIV Bl The value for the cutoff A has been fixed 
to 600 MeV, and d R = 1.75 and g 4V = 0.793. Notation is the 
same as in Table rVTHl 

on realistic interactions and weak currents consisting 
of vector and axial-vector components with one- and 
many-body terms. Two different approaches have been 
used to derive these operators: the first one goes be- 
yond the impulse approximation, by including meson- 
exchange current contributions and terms arising from 
the excitation of A-isobar degrees of freedom. This ap- 
proach, labeled SNPA, has been widely and successfully 
used in studies of electroweak processes (see for instance 
Refs. [20, [2l|, IS^j S^l ) • The second approach, labeled 
EFT*, includes two-body contributions, beyond the im- 
pulse approximation, derived within a systematic %EFT 
expansion, up to N3LO. The only parameter in the SNPA 
nuclear weak current model is present in the axial cur- 
rent and is determined by fitting the experimental value 
for the triton half-life. In the case of the EFT* ap- 
proach, two LECs appear, one in the vector and one in 
the axial-vector component (note that the LEC appear- 
ing in the isoscalar electromagnetic contact term does 
not contribute here). They are fixed to reproduce, re- 
spectively, the A = 3 isovector magnetic moment and 
triton half-life. 

Our final results are summarized in Eqs. (15. ip 
and (15. 2|) . The very accurate experimental datum of 
Ref. [1J| for the total rate in muon capture on 3 He is 
very well reproduced. For the muon capture on deuteron, 
a precise measurement should become available in the 
near future @ . The dependence of the results on the in- 
put Hamiltonian model, or on the model for the nuclear 
transition operator, is weak, at less than 1 % level. This 
weak model dependence is a consequence of the proce- 
dure adopted to constrain the weak current. 

We conclude by noting that (i) within the EFT* ap- 
proach, the two-body vector current operators beyond 
the one-pion-exchange term give significant contributions 
to the total rate, especially for muon capture on 3 He. 



These terms have been included in the study of weak 
processes here for the first time. However, since they are 
proportional to the momentum transfer, their contribu- 
tion is expected to be negligible in the pp and hep reac- 
tions [22J. (ii) Some of the radiative corrections are ac- 
counted for in the present predictions for the muon rates, 
since the value adopted for the Fermi coupling constant 
Gy is that extracted from an analysis of superallowed 
/3-decays (44[. The size of these corrections (roughly 2.4 
% on G v ) is i n ag reement with that calculated indepen- 
dently in Ref. [7l( for the muon capture on hydrogen and 
helium. The additional radiative corrections, originating 
from vacuum polarization effects on the muon bound- 
state wave function, have also been estimated by the au- 
thors of Ref. [7l[ . In the case of the muon capture on 3 He, 
they increase the predicted rates by about 0.68%, leading 
to 1496 s" 1 for the AV18/UIX model, and to 1490 s" 1 
for the EFT*(N3LO/N2LO) model, and thus bringing 
them into closer agreement with experiment, (iii) The 
value for the induced pseudoscalar coupling used in the 
present study, gps — —8.28 at the four-momentum trans- 
fer relevant for muon capture on 3 He (q% = 0.954m^), is 
consistent with that predicted by PCAC and chiral per- 
turbation theory [72j . The agreement between the calcu- 
lated and experimental muon capture rates confirms the 
validity of these predictions, (iv) The agreement between 
our final results for the muon capture rate on 3 He and 
those of Ref. [13] confirms the tight limits obtained there 
on possible contributions from second class currents. 

Finally, we remark that it would be interesting to ex- 
tend these calculations to the processes fi~ + 3 He — >• 
n + d + and fj,~ + 4 Hc — > n + 3 H + f M , for which 
experimental data are also available [ii, [74| . 
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